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FOREWORD 


The  MITRE  Corporation  is  concerned  with  the  survivability  of  the 
Air  Force  Command  and  Control  Systems.  It  conducts  studies  in  this 
general  area  in  order  to  determine  the  levels  at  which  various  systems 
components  fail  and  investigates  various  alleviating  measures  which  may 
be  employed  to  raise  the  levels  of  survivability. 

One  phase  of  this  work  is  concerned  with  the  behavior  of  deep  under¬ 
ground  hard  command  posts  excavated  in  soil  and  rock  when  subjected  to 
nuclear  attack. 

Among  the  many  problems  involved  in  this  area,  one  of  the  most 
important  and,  perhaps,  least  understood  is  the  phenomena  associated  with 
the  transmission  of  ground  shock  to  the  underground  command  post.  It  is 
known  that  shock  loads  such  as  those  produced  by  nuclear  weapons  will  be 
transmitted  through  the  ground  by  means  of  stress  waves;  however,  because 
of  the  complexity  involved,  only  skeletal  information  is  presently  avail¬ 
able  that  can  be  applied  to  the  actual  design  of  the  underground  command 
post.  It  is  important  that  basic  research  in  this  field  be  accelerated 
so  that  appropriate  criteria  can  be  established  for  the  design  of  installa¬ 
tions  of  interest  to  the  Ccmmand-Control  Development  Division. 

At  MITRE  programs  have  been  initiated  to  study  stress  wave  propagation 
in  various  media  because  we  are  interested  in  many  given  geographic  locations, 
and  since  earth  materials  which  vary  widely  from  site  to  site  exhibit  different 
properties  and  characteristics,  many  analytical  models  are  necessary  to  pre¬ 
dict  the  response  behavior  of  the  geologic  materials  at  a  variety  of  sites. 

Some  of  these  models  are:  linear -elastic,  non-linear  elastic,  elastic-plastic, 
elastic -locking,  visco-elastic,  visco-plastic,  etc. 

Information  contained  herein  is  concerned  with  stress  wave  propagation 
phenomena  in  visco-elastic  media.  The  numerical  results  given  are  for  the 
purpose  of  demonstrating  visco-elastic  effects  and  should  not  be  construed 
as  representing  the  actual  behavior  of  a  given  rock  or  soil  medium.  This  is 
due  to  a  lack  of  information  on  the  physical  visco-elastic  constants 
associated  with  various  types  of  rock  and  soil. 

This  report  is  a  part  of  a  series  of  studies  currently  being  carried  out 
by  MITRE  and  Raul  Weidlinger,  Consulting  Engineer,  New  York  City. 
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INTRODUCTION 


This  report  Is  one  of  a  series  of  technical  discussions  and  papers 
concerned  vith  the  theory  of  vave  propagation  in  solids  with  special 
applications  to  ground  shock  phenomena.  It  presents  theoretical  results, 
on  the  free  field  effects  due  to  progressing  pressure  loadings  on  the 
surface  of  a  semi -infinite  linearly  viscoelastic  (standard  solid)  half- 
space. 

The  ultimate  purpose  of  this  group  of  papers  is  to  arrive  at  conclusions 
for  the  free  field  effects  due  to  progressing  surface  pressure  loadings  on 
actual  media  such  as  rock  or  soil  vith  coup  lex  properties  vhlch  are  difficult 
to  analyse.  The  problem  of  rocklike  media  has  been  approached  by  considering 
a  succession  of  materials  having  gradually  more  complex  properties  such  as 
1)  an  acoustic  lnvlsdd  fluid,  2)  a  linear  elastic  solid,  3)  a  non-line arly 
elastic  solid  £(l)-(3)  -  See  Reference  [1]J  and  a  linearly  viscoelastic  solid 
In  the  present  report.  It  Is  felt  that  In  this  manner,  certain  conclusions 
vhlch  can  be  drum  for  cases  of  simple  properties  can  be  extrapolated  for 
more  complex  properties  by  qualitative  reasoning. 

The  paper  represents  a  step  In  the  investigation  of  free  field  phenomena 
In  viscoelastic  materials  but  considerable  future  vork  remains  to  be  performed 
in  this  field.  It  presents  the  plane  strain  solution  for  the  stress  distribution 
produced  by  the  uniform  motion  of  the  pressure  vave  on  the  surface  of  the  medium 
The  analysis  is  based  on  the  assumption  that  a  steady  state  exists  vith  respect 
to  a  coordinate  system  attached  to  the  moving  load  and  that  the  velocity  of  the 
moving  load  is  greater  than  that  of  lrrotatlonal  and  equivoluminal  vaves  in  the 
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medium  (superselsmie  cam),  further  Investigations  in  which  the  free  field 
streeeea  in  the  trans-seisaic  and  eubeeleaic  regions  are  determined  and  In 
which  the  velocity  of  the  surface  pressure  varies  as  a  function  of  time 
will  be  required  for  a  more  complete  picture  of  the  phenomena  In  question* 

Attention  stust  also  be  paid  to  the  relation  of  the  mathematical  model 
to  the  expected  physical  behavior  of  the  material,  for  those  rocks  which 
nay  be  expected  to  act  vlscoelastlcally,  much  work  remains  to  be  done  on 
the  experimental  determination  of  the  appropriate  viscoelastic  constants 
which  are  required  as  Input  parameters  In  the  analysis]  such  Information 
is  not  available  at  the  present  time. 

Numerical  results  are  presented  for  a  hypothetical  viscoelastic  material 
and  the  free  field  stresses  are  evaluated.  They  are  compared  with  the  cor¬ 
responding  stresses  in  the  material  in  Its  relaxed  and  unrelaxed  elastic 
state.  The  results  are  Illustrative  only  and  should  not  be  taken  as. 
applying  to  specific  real  materials  unless  experimental  Investigations 
show  the  coincidence  of  the  viscoelastic  parameters  for  the  modal  and  the 


real  material 
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*•  g  2  -  - 

At  Viscoelastic  operator:  A  -  (5^  -  1)  ♦  • 

X  E 

:  Nondlmentlon&l  space  coordinates:  I  -  ^  *  H  ■  • 

i  Ratio  of  propagation  velocity  of  high  frequency  equivoluainal 

waves  to  high  frequency  Irrotational  waves.  X  •  ^  gjj  • 

X  (X  ):  Unrelaxed  (relaxed)  Lam£  constant. 

fu  ^  «v 

X  :  Viscoelastic  Lame  operator.  X  »  X  -  ^  p  . 

41 

p  ()i  ):  Unrelaxed  (relaxed)  shear  modulus, 

p  :  Viscoelastic  shear  operator.  £  ■  p 

v:  Poisson's  ratio  for  the  unrelaxed  material. 

p:  Density  of  the  material. 

o^j  :  Cartesian  components  of  the  stress  tensor. 

(„  ).  Portion  of  stress  components  due  to  the  irrotational 

(equivoluainal)  potential. 


loalor  (or  lrrotationol)  potential. 
Vector  (or  equieoluninal )  potential. 

The  y  oavonent  of  • 
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latroductlon. 

This  paper  coaaidera  the  plane  strain  problem  of  finding  the  stress 
distribution  produced  by  the  uniform  motion  of  a  step  pressure  on  the 
surface  of  a  linear  viscoelastic  (standard* solid)  half-space  (Fig.  l). 

Ihs  speed  of  the  load  is  assumed  to  be  greater  than  the  velocities  of 
plane  lrrotational  and  equivoluminal  vaves  in  the  medium  (super seismic  case), 
eliminating  the  effect  of  initial  conditions,  it  is  assumed  that  the  load 
is  moving  in  from  x  •  »,  such  that  a  steady  state  exists  with  respect  to 
a  coordinate  system  attached  to  the  moving  load.  The  equivalent  problem 
for  the  elastic  case  has  been  treated  in  Ref.  [1). 

The  solution  for  arbitrary  load  distribution  can  be  obtained  by 
superposition  from  the  solution  for  the  step  pressure.  Although  not 
considered  here,  the  case  of  tangential  loads  on  the  surface  of  the  half- 
space  could  be  treated  in  the  same  manner. 

Formulation. 

I*t  (x,  y,  i)  denote  fixed,  rectangular,  Cartesian,  space  coordinates, 
vhile  (x,  y,  s)  denote  space  coordinates  attached  to  the  front  of  the 
moving  load  pQ  H(x)  (Fig.  1),  where  pQ  is  the  pressure  intensity  and  H  is 
the  unit  step  function.  The  uniformly  distributed  pressure  moves  over 
the  surface  of  the  half-space  in  the  negative  x  direction  with  a  speed  V. 

At  the  time  t  •  0,  the  two  coordinate  systems  are  taken  to  be  in  coincidence 
so  that 


x  •  x  ♦  Vt 


#  y  -  y  > 


i  •  i 


(i) 
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The  stress -displacement  relations  in  a  homogeneous,  Isotropic  medium 
which  is  elastic  in  bulk  and. viscoelastic  in  shear  is 


•u 


x  V  *  ^ui,j  ♦  uj,i> 
i,J,k  -  x,y,i 


(2) 


with 


(3.) 


where  K  is  the  bulk  modulus  and  p  is  the  shear  operator?  Selecting  the 
standard- solid  model: 


E  -  u 


1  +  IS 


(3b) 


where  p  la  the  "unrelaxed "  and  mp  ■  p*  (0  <  a  <  1.)  is  the  relaxed  shear 
modulus,  while  T  is  a  relaxation  time  (see  Ref.  2). 


The  equations  of  motion  of  a  continuum  in  the  absence  of  body  forces 
are,  for  small  displacements, 


°U,J  "  pSi 

where  p  is  the  density.  Eqs.  (2)  and  (4)  can  be  combined  to  give  the 
differential  equations  on  the  displacements: 

(*  ♦  S>uj,ji  ♦  5  ui,jj  •  p«i 


<^) 


(5) 


By  means  of  the  Helmholts  resolution,  the  displacement  vector  can  be 
separated  into  an  irrotational  and  equivoluminal  part 

'  tw-0  (6) 

where  is  the  alternating  tensor,  and  where  the  irrotational  part  comes 

from  the  scalar  potential  }  and  the  equivoluminal  part  from  the  vector 

,  / 

potential 

*  The  usual  convention  of  summation  over  repeated  subscripts  is  adopted. 
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Equation  (5)  will  be  satisfied  If  j  and  ^  satisfy  the  following  equations} 


(K  +  2m)  i  ,J(j  -  p  1 

(7*) 

•  • 

S*i,jj  -  p*i 

(7b) 

where 


^  d  +  Tc 

K  ♦  2m  -  (X  +  2M)  ( - £) 

11+  *3fc 

X-K-Sjj 


n  *  1  -  ^^(1-m) 

V  .  — t _ _  4=gv 

-*■  X  +  2m  2(1- v) 


(7C) 


v  being  Poisson's  ratio  for  the  unrelaxed  material.  Rote  that  for 


-1 <  y  <  i  »  0  <  n  <  1. 


For  plane  strain  in  the  (x,  £)  plane,  a  suitable  torn  for  the  potentials 

is 

$-  §<x,  l,  t)  ,  fj  -  -  0  ,  i-  -  ^(i,  l,  t)  (0) 

Since  the  solution  is  a  steady  state  solution  in  the  (x,  s)  coordinate 
system,  transforming  to  this  system  leads  to  equation  Independent  of  t. 

Thus,  by  the  use  of  Eqs.  (1)  and  (0),  Eqs.  (7a)  and  (7b)  give: 


(9»>) 


■HdNMWfcW*  •• 


where 


C?  .iiJLiH 

L  p 

o  *  • 


4  • 

p  L 


<$.M 

T  p 


#2  * 

c!  .  n- 

T  p 


mc£ 

T 


(9c) 


I  *  W  1  n  "  VT  J 

Hie  quantity  (C^)  represents  the  velocity  of  prop  station  of  high  (low) 
frequency,  plane,  lrrotatlonal  vaves  In  the  standard- solid  medium,  while 
CT  (C*)  represents  the  propagation  velocity  for  high  (low)  frequency,  plane, 
equlvoluminal  waves. 


The  boundary  conditions  to  be  satisfied  by  £  and  i  are  determined 
from  the  traction  conditions  at  the  surface  of  tbs  half  space,  q  ■  0  t 


»„(l*°)  ■  0  (10a) 

*„(l,0)  -  -  poH(|)  (lCb) 


Utilising  Iqs.  (1),  (2),  (6)  and  (9),  Kqs.  (10)  can  be  written  In  terms  of 
Jand  -J  .  At  i)  •  0  1 


(Ua) 

(Ub) 

(11c) 


-13- 

Tfce  formulation  of  the  problem  le  not  yet  complete,  for  there  le 
•till  a  further  condition  to  be  specified  vhlch  arises  from  physical 
considerations.  Because  the  load  moveii  with  a  velocity  which  le  greater 
than  the  propagation  velocities  in  the  standard* solid  medium,  no  disturbances 
can  ever  get  ahead  of  the  load.  Thus  for  |  <  0,  the  medium  is  undisturbed 
(see  fig.  (l))«  Actually,  since  the  propagation  velocities  of  lrrotatlonal 
waves  have  a  definite  maximum  value  (C^),  it  is  seen  physically  that  there 
will  be  a  straight  line  of  demarcation  between  the  disturbed  and  undisturbed 
portion  of  the  medium  *  the  Mach  line  (or  wave  front)  for  lrrotatlonal 
disturbances,  indicated  on  Fig.  (1)  by  the  line  01.  lrrotatlonal  disturbances 
can  exist  only  behind  this  line.  Similarly,  equlvolumlnal  disturbances  can 
exist  only  behind  the  equlvolumlnal  wave  front  or  Mach  line,  OE  In  Fig.  (l). 
Thus  for  |  <  0  quiescent  conditions  exist,  and  the  solution  will  have 
non- vanishing  values  only  for  |  >  0,  permitting  the  use  of  a  Laplace  transform 
in  (• 

The  formulation  of  the  problem  is  now  complete .  Equations  (9)  nust.be 
solved  subject  to  the  conditions  at  q  ■  0,  given  by  Bqs.  (11),  and  to  the 
condition  of  quiescence  for  q  <  0.  Having  obtained  $  and  $ ,  the  displacements 

a 

and  stresses  may  be  computed  from  Eqs.  (6)  and  (2),  respectively,  utilising 
Eqs.  (l)  and  (8). 

Formal  Solution. 

The  solution  to  the  problem  can  be  written  formally  as  a  complex 
integral  by  the  use  of  the  Laplace  transform,  defined  by 

•0 

*(p»tj)  ■  J  di  (12) 


with  the  inversion 


"'•’'IT  -'V* - 


-lU- 


7+i— 


f(l#n)  •  2?r  y  *( p.nV*  <Jp 


(13) 
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Applying  this  transform,  and  utilising  the  condition  of  quiescence 
for  |  <  0,  Eqs.  (9)  yield 


d2f  _2-  2 


-  p\2  i  -  0 


d2$  2-  2  j 

Tjr  -  p  “t  ’  •  0 

OT| 


where 


■l-  "L 


L  P+n 


r- 


Re  a.  >  0 


®T  “  “t  LnF* 


« 

p«,<Vv2 


>  0 


»L  -  <4  ■ 1)1/2  I  “r  -  *  Dl/2 


-£-<*£  .l)1/zi^(<2.l) 


KL  -  JL  .  V  ,  »»*  V 

^  C,  '  ^  “  7  ^"7 

L  CL  CT 


Equations  (l4a,  b)  have  the  solutions 


1/2 


I . 

,T  -P“mn  PV1 
T  -  Be  ♦  B'e  1 


(1W) 

a1*) 


(14c) 


(15*0 


(15b) 


where  A,  A',  B,  B'  are  arbitrary  constants.  The  condition  of  quiescence 
for  |  <  0  will  be  satisfied  if 


A*  -  B'  -  0 


(15c) 


. - 
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The  constants  A  and  B  are  determined  from  the  transforms  of  the  boundary 
conditions,  Eqs.  (ll): 


2  ^  -  pfij2  -  1)  i  ■  0 

^  -  1)1  ♦  2  I-  *  (VT)2  (~ 


These  equations  give  for  A  and  Bt 


rp0tvT)ai  r 

5+1  1  hj  il 

L  mA  Jl 

■p3(p+ar  •*  ■* 

B  -  2 


where 


A  -  (5j,2  -  1  f+  ^  ijij 


(16a) 

(16b) 


(17*) 


(17b) 


(17c) 


Thus  T  and  ^  are  determined,  and  the  transforms  of  the  stress  components 
can  now  be  found. 


8treaa  Components  at  the  Wave  fronts. 

From  Eqs.  (2)  and  (6)  It  is  seen  that  the  stress  consonants  any  be 
written  In  two  portions 


”u  ■  "u,  *  "u* 


(16) 


where  results  from  the  lrrotatlonal  potential  § 
equlvolualnal  potential  'f'. 


otJt  from  the 
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Ut ill ting  the  results  of  ttaa  preceding  sections,  the  transforms  of  thee* 
quantities  are  t 


.  v]  [■,*  -  4*' 

(1?*) 

k-  Wh‘4’1' 

(19b) 

fc]fc«r2 

- 

(19c) 

-  »}-** 

« 

(19d) 

2  - 
-PVI 

s 

(19«) 

• 

(19f) 

-  -  fc]  IV  -  -  *]  [*«r2  • 

(19«) 

0 

(19h) 

Note  the  formal  correspondence  between  these  results  and  the  results  for 
the  elastic  body.  If  the  Laplace  transform  (In  the  variable  |  ■  of  the 
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elastic  solution  (Rsf .  (1)  -  Eqs.  (51-53))  i*  then  ^  and 

in  the  elastic  solution  are  replaced  by  m^  end  ,  Eqs.  (19)  obtained. 

Because  of  the  complexity  of  the  inverted  forms  of  and  *  It 

is  useful  to  obtain  closed  form  expressions  for  these  quantities  and  their 
first  derivatives  at  the  wave  fronts.  These  expressions  are  obtained  by 

PV»  - 

considering  the  asymptotic  expansions  of  the  quantities  e  and 

J>«y\  ^  A>  ghown  la  the  Appendix,  the  first  two  coefficients  of  the 

expansion  in  powers  of  ±  of  e  c1J(f  represent  the  values  of  and 

•it  .  Just  behind  the  irrotational  wave  front.  A  similar  situation 
d|  ijp  " 

exists  for  the  equlvoluminal  portion  of  the  stress  consonants  at  the 
equivoluminal  wave  front. 

Isssediately  behind  the  irrotational  wave  front  (i.e.  at  |  »  s^n+)» 
these  values  are: 


a-~*  -  -  <£)  (i^  -  (-4  ' 

*0 


(20a) 


(20b) 


(20c) 


(20d) 
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a  .  i§(4  -  1)  (1  -  -)  ♦  ^  4  (1  *  n)  +  4  U  ’  *} 
b  *  [b5^1  [^(1+3n)  ■  bn]  ® 


Immediately  behind  the  equivoluminal  wave  front  (i.e.  at  1  -  ) 

the  vaiua  and  first  derivatives  of  arei 


(201) 


gxx+  Vt  #-Pn  (21a) 

*o 


(21b) 
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The  quantities  and  are  the  equivalent  Jumps  In  an  elaatle 
medium  having  elastic  moduli  equal  to  the  unrelaxed  moduli  of  the 
standard  solid. 

1 

Stress  Components  Par  Behind  the  Wave  fronts. 

To  obtain  the  value  of  the  stress  coeponents  far  behind  the  vave 
front,  the  final  value  theorem  for  the  Tnplace  transform  is  utilised  t 

lln  p5u.lU  oi(} 

H* 

The  values  obtained  are  (for  |  -*  »)  i 


(22*) 

ia.Si 

Po  ^ 

(22b) 

°xt»  ^®L  ,  . 

-  Hf  («*  •  i) 

¥o  b 

(22o  ) 

0  .  m  -  0 

xsf  xsp 

(22d) 
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Po  ’  b* 
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^•-(7)(-v^)(<2-0  <“«> 

0^  •  0  (22h) 

where 

6*  -  (*£-  1)  +  (221) 

These  values  correspond  to  the  values  of  the  stress  components  which  would 
occur  in  an  elastic  medium  having  elastic  moduli  equal  to  those  of  the 
totally  relaxed  standard-solid*  This  fact  might  have  been  foreseen  from 
purely  physical  reasoning. 


iroximate  Expressions  for  the  Stress  Components. 


By  utilising  the  egressions  for  the  Jumps  and  first  derivatives  of 
the  stresses  behind  the  wave  fronts,  and  for  the  values  of  the  stresses 
as  |  •*  •»,  approximate  expressions  for  the  stress  components  for  all 


values  of  |  may  be  constructed,  for  Instance,  having  the  value  of  o. .  and 

ijq> 

at  the  irrotational  wave  front,  and  of  as  |  -»  «,  a  suitable 
curve  may  be  Interpolated  to  connect  the  value  of  o. ,  (with  proper  slope) 
at  the  wave  front  and  the  asymptotic  value  of  at  |  An  exponential 

type  of  response  is  typical  of  linear  viscoelastic  phenomena,  and  such  a 


curve  if  therefore  suggested.  This  suggestion  is  supported  in  Ref.  [3] 


where  it  is  shown  that  such  an  approximation  is  succesful  in  an  analogous 
problem  (wave  propagation  in  a  rod  of  standard-solid  behavior). 
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Appendlx 


To  show  that  the  values  and  derivatives  of  immediately 
behind  the  lrrotational  wave  front  can  be  determined  from  the 

pawn  _ 

coefficients  of  the  asymptotic  expansion  of  e 
Heaviside's  series  expansion  method  is  employed  (Ref.  (4]). 

P«V n 

Suppose  e  5  >  (p,T|)  can  be  expanded  in  the  series  j 

yp 


pawn 


m 


(n)  p 


-n 


Inverting  both  sides 


«iJ(p<i  ♦  -  H(t) 


f  *.<i)  i”'1 
L  u-i« 


where  the  right  hand  side  has  been  Inverted  term  by  term  (Justified  in 
Ref.  (4} ),  and  the  left  band  side  has  been  inverted  by  use  of  the  shift 
theorem  with  the  knowledge  that  *  0  for  |  <  m^V  Thus 


°ij*  (M) 


h(i  -  m^)^ 


An(n)  (t-Mj.n)1 


n-l 


n»l 


But  the  infinite  sum  represents  the  Taylor  series  expansion  in  (,  about 
|  ■  m^Tj+,  of  (continued  analytically  at  |  ■  m^t)+). 

Hence 


d|n  W 
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Problema  involving  stress  wave  propagation  in  viscoelastic  materials 
can  generally  be  treated  by  Integral  transform  methods,  but  the  complete 
numerical  evaluation  of  the  inversion  Integral  becomes  in  many  cases 
very  complex  and  lengthy.  It  is  usually  possible,  however,  to  use 
asymptotic  methods  to  Obtain  the  value  of  the  stress  c(tQ)  and  its  time 
derivative  e(tQ)  at  the  time  of  arrival,  tQ,  of  the  stress  wave  at  a 
particular  point  in  the  body}  in  many  cases,  these  initial  values  may  be 
obtained  in  closed  form.  In  addition,  the  long  term  solution  for  the 
stress,  o(«),  can  also  be  determined  in  a  simple  manner,  either  from  the 
physical  situation  or  from  an  asymptotic  evaluation  of  the  inversion 
Integral. 


Once  the  values  of  o( tQ),  d(tQ)  and  c(»)  have  been  found.  It  is 
possible  to  Interpolate  a  curve  for  the  stress  o.  Bleich  and  ** 
have  presented  an  exponential  interpolation  in  Reference  [1] .  Their 
interpolation  is  based  on  the  reasoning  that  since  an  exponential  type 
of  response  is  typical  for  linear  viscoelastic  media,  an  exponential 
interpolation  is  appropriate.  Figures  (la,  b)  show  typical  stresses  a 
(at  some  point  in  the  body)  to  which  the  interpolation  may  be  applied 
as  a  function  of  an  actual  or  non-dimensional  time  t.  The  complete 
history  of  the  stress  is  given  by  E4.  (1)  of  Reference  [1): 
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vhare  t  ,  o  .  a  and  a  are  evaluated  first  for  the  problem  under 
o  o  o  • 

consideration . 


In  the  application  of  the  above  interpolation  to  practical  problem# 
it  is  also  possible  to  encounter  cases  in  which  (am  -  o^)  is  opposite  in 
sign  to  Cq  •  For  such  cases#  the  exponential  interpolation  of  Eq.  (1) 
is  obviously  invalid  since  the  exponential  would  be  raised  to  a  positive 
power  in  time*  Figures  (2a#  b)  show  typical  stress-time  relations  of 
this  type.  An  interpolation  in  the  form  of  the  product  of  a  linear 
polynomial  in  t  multiplied  by  an  exponential  in  t  can  be  used  in  such 
a  case.  The  complete  history  of  the  streos  is  then  given  byt 


t  <  t 


-  [(°-  -  •o)(1  +  IT2)  *  Ao<t't0>]  •  ^ 


t  >t> 
—  o. 


1(2) 


In  the  case  of  Eq.  (1)#  the  interpolation  was  uniquely  determined  by  the 
three  values  oQ#  aQ  and  .  In  the  present  interpolation,  Eq.  (2)# 
these  three  quantities  and  an  additional  constant  T  are  required. 

The  value  of  T  cannot  be  determined  unless  an  additional  asymptotic 
value  is  obtained  from  the  inversion  Integral  of  the  problem;  this  in 
general  would  be  very  difficult  to  do.  It  suggests  itself#  however#  to 
chose  the  decay  time  of  the  viscoelastic  medium  as  a  suitable  value  of 
the  constant  T. 


The  interpolations  of  Eqs.  (l)  and  (2)  may  be  utilised  to  obtain 
etreoses  which  are  produced  by  tine  varying  pressure  loadings  P(t)  acting 
on  the  viscoelastic  body#  once  the  stresses  which  are  produced  by  a  unit 
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step  pressure  loading  are  evaluated.  Let  the  quantities  e0«  ffQ  and  em 
be  the  stresaes  at  a  point  in  the  viscoelastic  body  which  are  produced 
by  a  unit  step  pressure  loading.  The  corresponding  stress-tine  history 
c(t)  which  is  produced  by  the  loading  P(t)  is  obtained  through  the  use 
of  Duhamel's  Integral  and  Eq.  (1)  or  Eq.  (2). 


for  an  applied  pressure  of  the  form 

at  St 

m  mb  m 

5  +  *«  p]  (3) 

the  stress  e(t)  produced  at  a  point  in  the  body  becomes,  using  Eq.  (1): 


o(x)  ■  0  t  <  0 


where  t  -  t-tQ  is  measured  fron  the  tine  of  arrival,  %Q,  at  the  wave  at 
the  point  and 


(5). 


Repeating  the  procedure  with  Eq.  (2),  the  stress  at  a  point  in  the 
body  due  to  the  tine  decaying  pressure  loading  of  Eq.  (3)  become  si 
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where 


A  "  °m  *  «o 


and 


(*) 

0>) 


(7) 


The  Interpolations  considered  should  prove  convenient  in  different 
problems  of  wave  propagation  in  viscoelastic  media.  It  can  be  applied 
in  two  -  or  three  dimensional  problems  where  two  distinct  signals,  similar 
to  P  and  8  waves  in  elastic  media  are  received  at  a  point.  In  such  cases, 
each  component  of  the  signal, may  be  approximated  by  an  expression  similar 
to  those  given  by  Eq.  (1),  (2),  (4)  or  (6).  The  applications  of  these 
Interpolations  to  the  stresses  produced  by  a  progressing  surface  pressure 
on  the  surface  of  a  viscoelastic  half-space  [Section  I  of  this  report] 
are  given  in  Section  XII. 
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«YALUATIO«  car  THE  BBUE88E8  IW  A  VISCOELASTIC  HALT-SPAC1  WHICH  ACT  FBOPUCIP 
BT  A  HPBBMgB  TDCS-IBCATIBO  OTTAB  FCT86UCT. 

1.  Introduction. 

The  stresses  which  arc  produced  in  a  linear  viscoelastlo  (standard 
•olid)  half- apace  by  the  uni  fora  motion  of  a  step  pressure  on  its  surface 
bars  been  studied  is  Section  X  for  the  case  in  which  the  velocity  of  the 
traveling  surface  pressure  is  greater  than  the  velocities  of  plane  lrro- 
tational  and  equivoluainal  waves  in  the  medlua  (superseisale  case). 

Specific  results  were  presented  which  allow  the  determination  of  the 
quantities  cq,  iQ  and  am  corresponding  to  both  the  lrrotatlonal  and  the 
equlvolumlnal  waves  produced  in  the  nediua  by  the  unit  step  pressure  on 
the  surface.  The  quantities  aQ,  aQ  and  oa  can  be  used  in  conjunction  with  . 
the  interpolations  which  are  presented  in  Section  II  to  determine  the 
stress  history  at  a  point  in  the  medium  produced  by  a  time-decaying 
surface  pressure  [Bqs.  (h)  or  (6),  Section  II). 

The  present  Section  presents  numerical  results  for  the  stress  cosponents 
9xx*  °»  1104  °»  at  various  points  in  the  viscoelastic  medium  [rig.  (l)] 
which  are  produced  by  the  time-decaying  surface  pressure  shown  in  Fig.  (2). 
(See  Eq.  (3)»  8ectlon  IX).  For  comparison  purposes,  the  corresponding 
stresses  which  would  be  produced  in  an  elastic  medium  having  a  shear  modulus 
equal  to  the  unrelaxed  shear  modulus  of  tbs  viscoelastic  medium  are  also 
shewn  in  aach  ease. 


'.'♦Nv'asekilrMk..  »'  • 
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The  stresses  produced  by  the  irrotational  end  the  equivoluainal  waves 
In  the  viscoelastic  media  ere  first  shown  separately.  In  addition,  the 
total  stresses  at  a  point  are  evaluated  by  .  -  superimposing  the  irrota¬ 
tional  and  the  eq,uivoluminal  stresses  with  an  appropriate  tine  delay. 

Due  to  the  absence  of  reliable  information  on  the  appropriate  visco¬ 
elastic  constants  for  reed  materials,  a  set  of  viscoelastic  constants  has 
been  chosen  for  an  illustrative  example  only.  The  results  should  not  be 
taken  as  applying  to  a  specific  real  material  unless  the  coincidence  of 
the  viscoelastic  parameters  for  the  real  material  and  those  used  in  the 
illustrative  example  can  be  shown.  The  numerical  values  of  the  material 
and  the  load  constants  which  were  used  in  the  computations  are  tabulated 
below. 

(a)  Material  Constanta  -  Viscoelastic  Medium. 

|i  ■  Uhrelaxed  shear  modulus  •  1.200  (10^)  lb/in^. 

M*  "  Relaxed  shear  modulus  ■  p/2  ■  0.600  (10^)  lb/in2, 
p  ■  Mass  density  of  medium  ■  $.166  lb, see2/ft\ 

(v  •  167  lb/ffc3). 

v  •  Poisson's  ratio  of  unrelaxed  body  •  0.2$. 

T  a  Relaxation  time  for  viscoelastic  medium  ■  1$  ms. 

(b)  Material  Constants  -  Linearly  Elastic  Media. 

p  ■  Mass  density  of  medium  ■  $.186  lb.sec^/ft** 
v  *  Poisson's  ratio  ■  0.2$. 

0^  ■  Telocity  of  P  waves  ■  10,000  ft/sec  [unrelaxed  p). 

Cg  ■  Velocity  of  8  waves  ■  6,000  ft/sec  [unrelaxed  p]. 

p  ■  Shear  modulus  -  p  (unrelaxed)  ■  1.2  (10^)  lb/in2. 

or 

p  ■  Shear  modulus  a  p  (relaxed)  ■  0.6  (10^)  lb/in2. 
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Tha  alastie  medium  for  which  |i  (unrelated)  &■  used  will  ho  roforrod 
to  u  Klastle-tibrelaxtdj  that  for  which  p  (relaxed)  it  ueed  will  ho 
roforrod  to  at  Xlastle-Relaxed. 

(c)  Conotontt  for  Burfaco  Pretture  Distribution  [Fig.  (2)1. 

PQ  >  Peak  pressure  ■  8000  Vb/lt?. 

V  ■  Uniform  Telocity  of  tr art line  wave  «  12,000  ft/too. 

Relaxation  time  of  ourfaeo  proaturo  ■  30  at. 


2.  ttsmerlcal  Results. 

Conaidor  the  geometry  ahcnm  in  Pig.  (1)  in  which  tha  point  A  io 
located  at  a  depth  t  ■  300  ft  and  X  and  I  roproaont  the  irrotational  and 
the  equlvolumlnal  wave  front  a  reepectlvely  which  are  produced  hy  the 
aurface  proaturo  thorn  in  Pig.  (2). 

Figures  ( 3 )- ( 5 )  ahow  the  atreaaoa  and  ou  which  are  produced 

at  the  point  A  (1  ■  300  ft)  by  the  irrotational  wave  in  the  wiacoelaatie 
■odium,  the  value  of  tha  abscissa,  t  -  0,  roproaenta  tha  arrival  time 
of  the  wave  at  tha  point.  Tha  reeulta  wore  obtained  ueing  an  interpolation 
of  tha  type  ahown  in  Eq.  (4),  Section  ZXl  and  the  input  paranetara  «0,  am 
and  oQ  vara  cooputed  from  Eq.  (20)  and  Eq.  (22)  of  Section  X*  Xn  each 
cate,  tha  corresponding  stress  in  an  elastic  body  whoso  shear  modulus  la 
equal  to  tha  shear  modulus  of  the  unrelaxed  viscoelastic  body  la  alao  shown. 

Xt  it  aeen  that  the  high  peek  stress  which  la  predicted  hy  linear  ilaatlc  theory 
la  considerably  attenuated  by  the  viscoelastic  medium.  however,  at  later 
times,  the  atresaes  in  the  viscoelastic  body  may  he  higher  than  those  in  the 
linearly  elastic  body. 


The  corresponding  stresses  are  shown  in  Figs.  (5)-(6)  for  ths  csss 
of  ths  equlvolumlnal  wavs.  The  value  of  the  abscissa  t  •  0  represents 
ths  arrival  tins  of  ths  equlvolumlnal  ware  at  the  point  A.  Interpolations 
of  the  type  given  by  Kq.  (A)  were  used  in  evaluating  and  o%x. 

In  all  cases,  it  is  seen  that  the  peak  stresses  which  are  predicted  by 
linear  elastic  theory  are  significantly  attenuated  by  the  viscoelastic  medium. 
In  addition,  the  viscoelastic  stresses  are  lower  in  magnitude  than  the 
corresponding  elastic  stresses  over  the  entire  time  history. 

The  complete  stress-time  histories  which  are  produced  at  the  point  A 

can  be  evaluated  by  superimposing  the  stresses  shown  in  the  previous 

figures  with  an  appropriate  delay  time.  Figures  (7 )-(9)  show  the  total 

stresses  o _ ,  o__  and  o _ for  both  the  viscoelastic  and  the  elastic 

xx  ss  ax 

[unrelaxed]  media  where  the  time  v  ■  0  represents  the  arrival  time  of  the 
irrotatlonal  wave  I  at  the  point  in  question.  The  peak  stresses  which 
occur  at  the  significant  times,  l.e.  the  arrival  times  of  the  Irrotatlonal 
sad  equlvolumlnal  waves  at  the  point,  are  considerably  attenuated  by  the 
viscoelastic  medium.  Intermediate  values  of  the  stresses  mey  be  larger 
then  those  for  the  elastic  [unrelaxed]  medium. 

It  should  be  noted  that  the  significant  attenuations  of  peak  stresses 
by  the  viscoelastic  medium  are  purely  a  function  of  the  depth  a  and  do 
not  occur  for  relatively  shallow  depths.  To  illustrate  this.  Figs.  (lO)-(ll) 
show  the  magnitude  (not  the  direction)  of  the  principal  stresses  at  points 
located  at  depths  t  -  50  ft  and  t  ■  500  ft,  for  the  irrotatlonal  and  the 
equlvolumlnal  waves  respectively.  Interpolations  of  the  type  given  by 
Bq.  (4)  were  used  in  evaluating  e  for  the  Irrotatlonal  wave  while  those  of 
!q.  (6)  were  required  for  the  case  of  the  equlvolumlnal  wave.  For  the  stresses 
produced  by  the  Irrotatlonal  wave,  it  is  seen  that  the  attenuation  of  the  peak 
stresses  In  the  neighborhood  of  v  »  0  is  very  small  for  shallow  depths  but 
increases  significantly  with  depth. 
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For  tin  ratio*  v/C  encountered  In  this  example,  tha  vmlua*  of  tha 
P 

principal  lrrotational  alaatlo  atraaaa*  obtained  using  tha  unralaxad 
shaar  aodulus  arc  practically  aq.ua!  to  tha  eorraapooding  strasaaa  obtained 
by  using  tha  ralaxad  shasr  nodulus  and  consequently  only  *  single 
elastic  stress  curve  is  presented. 


Sea  Appendix  A ,  "Theoretical  Studies  in  Qround  Shock  Fhenouena",  by 
M.L.  Baron,  H.H.  Blaich,  P.  Valdlingnrj  The  NTERB  Corporatico, 
Report  SR-19,  October  i960. 


FIG. 2  SURFACE  PRESSURE-TIME  HISTORY 


STRESS  C£x  (IRROTAT  I ONAL  WAVE):  2  =  500  FT. 


STRESS  <r**  ( IRROTAT  I  ONAL  WAVE):  2*500FT. 


TOTAL  STRESS  Oxx  :  2=500FT. 


TOTAL  STRESS  (X,,  :  Z  =  500 FT. 


TOTAL  STRESS  OZ9  '■  Z  *  500  FT. 


FIG.  10  PRINCIPAL  STRESS-IRRO TATION A L  WAVE 

Z  =  50FT,  500FT. 
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